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A way to identify the would-be zero-modes of staggered lattice fermions away from the continuum 
limit is presented. Our approach also identifies the chiralities of these modes, and their index is 
seen to be determined by gauge field topology in accordance with the Index Theorem. The key idea 
is to consider the spectral flow of a certain hermitian version of the staggered Dirac operator. The 
staggered fermion index thus obtained can be used as a new way to assign the topological charge of 
lattice gauge fields. In a numerical study in U(l) backgrounds in 2 dimensions it is found to perform 
as well as the Wilson index while being computationally more efficient. It can also be expressed as 
the index of an overlap Dirac operator with a new staggered fermion kernel. 
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The computational efficiency and other attractive fea- 
tures of the staggered lattice fermion formulation made 
possible the first high-precision lattice simulations of 
Quantum Chromodynamics (QCD) [l|, and it continues 
to be widely used. However, this formulation has long 
been perceived as disadvantaged compared to Wilson 
fermions - the other traditional lattice fermion formu- 
lation - regarding topological aspects of QCD. This con- 
cerns in particular the Index Theorem 0] relating gauge 
field topological charge to fermionic zero-modes, which 
is needed to explain the large mass of the rf meson Q. 
A field-theoretic approach to the staggered fermion index 
was developed and used by Smit and Vink 0,0] ; however, 
it does not give an integer value for the index from the 
beginning and requires a renormalization that depends 
on the full ensemble of lattice gauge fields. In contrast, 
for Wilson fermions, an integer-valued index is obtained 
without any need for renormalization. This is because 
Wilson fermions have identifiable would-be zero-modes 
with definite chirality, provided the lattice gauge field is 
not too rough, from which the index can be defined. The 
Wilson fermion index 

0,01 , calculated from its equiva- 
lent description as the index of the overlap Dirac operator 
[j| , is currently widely used for issues such as calculation 
of the topological susceptibility [lfj and investigations of 
topological structure in the QCD vacuum - see, e.g 
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But for staggered fermions it has not been known how to 
identify the would-be zero modes and their chiralities in 
non-smooth gauge field backgrounds; they were only seen 
to emerge in simulations with improved actions that are 
sufficiently close to the continuum 

In this paper we show how the would-be chiral zero- 
modes of staggered fermions can be identified away from 
the continuum limit when the background gauge field is 
not too rough, thus determining an integer-valued index. 
A theoretical foundation for the Index Theorem for stag- 
gered fermions is established, placing them on the same 
footing as Wilson fermions in this regard. 



Our new approach in the staggered case parallels the 
spectral flow approach to the Index Theorem for contin- 
uum and Wilson lattice fermions 0, IH , which we begin 
by briefly reviewing in the following. Spacetime is taken 
to be a Euclidean box with periodic boundary conditions 
(the setting in which Lattice QCD simulations are per- 
formed) with even dimension d. Besides the physically 
relevant case d = 4 we will also consider d = 2 for illustra- 
tive purposes. From the hermitian Dirac gamma matri- 
ces {7^} M =i,...,d the chirality matrix 75 = -(i) d/2 7i • • • 7<j 
is denned. It has the properties 75 = 75, 7! = 1, 
{75 , 7m) = and {75 ,D}=0 where D = ^(d^ + A^) 
is the continuum massless Dirac operator on spinor fields 
coupled to a gauge field A. Consequently, the vectorspace 
of zero- modes of D (i.e. solutions to Dtp — 0) decomposes 
into ± chirality subspaces on which 75 = ±1. The index 
of D is the difference between the numbers n± of inde- 
pendent ± chirality zero-modes, and is fixed by topology: 
Gauge fields with smooth field tensor F fil/ (x) have an in- 
teger topological charge Q, and the Index Theorem in 
this setting states 



n+-n_ = {-l) d / 2 Q 



(1) 



The spectral flow perspective on the index arises by con- 
sidering the eigenvalues {A(to)} of the hermitian operator 



H(m) = 75(1? — m) 



(2) 



as a function of the parameter m. Note that a zero-mode 
tp of D with ± chirality is also an eigenmode of H(m) 
with eigenvalue A(to) = T m i crossing the origin with 
slope =Fl at m = 0. Furthermore, from the property 
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we see that these are the only eigenvalues of H(m) that 
cross the origin at any value of m. It follows that the 
spectral flow of H (to), defined as the net number of eigen- 
values A(to) of H(m) that cross the origin, counted with 
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sign ± depending on the slope of the crossing, comes 
entirely from eigenvalue crossings at to = and equals 
ri_ — n._|_, i.e. minus the index. 

In the lattice setting with Wilson fermions, the would- 
be zero-modes can be identified as the low-lying real 
eigenvalues of the Wilson-Dirac operator Dw 0. The 
spectral flow perspective 0, [I[ is based on the hermitian 
lattice analogue of ([2]): 

H w (m) = 75 (D w - to) (4) 

Regarding the spectral flows Hw(rn)ip(m) = X(m)ip(m), 
note that A(too) = Dwip( m o) = TOo^(too)- Thus 
eigenvalue crossings of Hw(m) are in one-to-one corre- 
spondence with real eigenvalues of D\y . Furthermore, af- 
ter normalizing eigcnmodes such that ip^tp = 1 one easily 
finds A'(to) = — -0(to)^75 , (/'(to) pj. Thus the sign of the 
slope of A(m) at a low-lying crossing value too is minus 
the chirality of the corresponding would-be zero-mode 
"0(too) for D\y. It follows that the index of Dw is minus 
the spectral flow of Hw{m) coming from the eigenvalues 
which cross the origin at low-lying values of to. Numeri- 
cal results illustrating this can be found, e.g., in Q- An 
illustration in the d — 2 case is given in Fig. 4 below. 

Turning now to staggered lattice fermions, where the 
lattice field x( x ) is scalar (rather than spinor) and de- 
scribes 2 d / 2 degenerate continuum fermion species (called 
quark tastes), the massless staggered Dirac operator is 

D st = V p (5) 

where V^ t is the usual lattice finite difference opera- 
tor coupled to the lattice gauge field and ij^xix) — 
(— f ) n i+---+ n w-i^-(x) where x = a(ni, . . . , rid) runs over 
the lattice sites. D st is anti-hermitian and therefore has 
purely imaginary spectrum. Consequently the identifica- 
tion of would-be zero-modes in the Wilson case does not 
carry over to the staggered case since it relied crucially 
on the role of non-zero real eigenvalues. For the same 
reason it is clear that an index of D st cannot be obtained 
from spectral flow of a staggered version of Hw{m). In 
fact the staggered analogue of (j4|), T 5 (D st — to), is not 
even hermitian. However, there is an alternative spectral 
flow approach which we now discuss, and which turns out 
to be perfectly suited to the staggered case. 

Return momentarily to the continuum setting and note 
that H(m) in ([2]) is not the only hermitian operator that 
can be used for the spectral flow perspective on the index. 
We could just as well use H(m) — iD — 77175. Its spectral 
flow is equal to minus the index of D just as before, since 
the previous argument, including the property (j3j), holds 
verbatim for this operator. But now the analogue in the 
staggered case, 

H st (m) = iD st - 771F5 (6) 

is also hermitian, so we can consider its spectral flow 
as well. Here T 5 is the analogue of 75 in the staggered 



formulation; it is hermitian and corresponds up to 0(a) 
discretization errors to 75 Si in the spin- flavor interpreta- 
tion [13]. Note that H st (0) — iD st , so the eigenvalues of 
D st are {— iA(0)} where {A(to)} are the eigenvalue flows 
of H st (m). This fact allows us to identify the would-be 
zero-modes of D st : as shown below, they are the eigen- 
modes with eigenvalues — iX = — iA(0) for which the as- 
sociated flow A(m) crosses zero at a low-lying value of to. 
Furthermore, the sign of the slope of the crossing is mi- 
nus the chirality of the would-be zero-mode, and hence 
the index is minus the spectral flow of H st (m) coming 
from the crossings at low-lying values of 777. 

To see this, consider first the situation in a smooth 
continuum-like gauge field background: The eigenvalues 
— iX of the would-be zero-modes x °f D st separate out 
from the rest of the spectrum; they are almost zero and 
have approximately definite chirality x^5X ~ if- F° r 
the corresponding eigenvalue flows A(to) of H st (m) we 
have A' (to) = — x( m ) ^5X( m ) j us t as m t ne Wilson case. 
Since x = x(0) it follows that A'(0) « +1, and since 
A = A(0) is almost zero it then follows that A (to) crosses 
the origin at a very small (in magnitude) value too of to. 
Whether too is positive or negative depends on the chi- 
rality and sign of A(0), but in either case the sign of the 
slope of the crossing is the same as that of A'(0), i.e. mi- 
nus the chirality. Under a roughening of the gauge field 
the location of a crossing will move, but it cannot disap- 
pear until it meets another crossing with opposite slope. 
Hence the would-be zero-modes remain identifiable with 
unchanged chiralities and index. 

To illustrate the identification of the staggered would- 
be zero-modes and their index we present results of a 
numerical study in U(l) backgrounds in 2 dimensions. 
Following [4} we start from specific smooth lattice gauge 
fields with topological charge Q and roughen them by 
multiplying the link variables by random phase factors: 
U fl (x) — > e^^^^U^x) with each r fl (x) randomly chosen 
in [— 7T, 71"]; the parameter S controls the roughness. Our 
numerical computations were all on the 12 x 12 lattice 
(as in Q) and done using ARPACK Q. The code was 
checked by reproducing the results in Table 1 of Q. 

Figure 1 shows low-magnitude eigenvalues of H st (m) 
versus m (horizontal axis) in a Q = 1 background of mod- 
erate roughness S = 0.33. The low-magnitude eigenval- 
ues of iD st are the eigenvalues at to = in the figure. 
There is no clear separation in their magnitudes. Nor 
is there a clear separation in the magnitudes of the chi- 
ralities as measured by x^sX'- numerical calculation of 
this quantity for the 3 pairs of eigenmodes with low- 
est magnitude eigenvalues gives —0.28, +0.17, +0.13. 
Nevertheless, the would-be zero-modes can now be iden- 
tified among these low-magnitude modes: as discussed 
above, they are precisely the eigenmodes whose eigenval- 
ues A = A(0) belong to eigenvalue flows A(to) which cross 
zero at a low-lying value of to. From the figure we see 
that there are two of these, both with positive slope, cor- 
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FIG. 1: Staggered spectral flow in a Q — l background. 



FIG. 3: Staggered spectral flow over a larger m range in same 
background as Fig. 1. 




FIG. 2: Staggered spectral flow in a 2 background. 

responding to negative chirality. Thus the index is —2. 
This is in accordance with the Index Theorem, since in 
the staggered fermion case ([1]) becomes 

index(C st ) = 2 d ' 2 (-\) d l 2 Q (7) 

As a second illustration, Figure 2 shows the eigenvalue 
flow in a Q = —2 background of the same roughness level 
6 = 0.33. From the eigenvalues at m = we see that 
iD st = H st (0) has one eigenvalue pair of much smaller 
magnitude than the others. Numerical calculation of 
X^r 5 x for these gives +0.34 compared to +0.21, +0.13 
for the eigenmode pairs with the next 2 lowest-magnitude 
eigenvalues. Hence we would naively expect the first pair 
to be the would-be zero-modes, giving index +2. How- 
ever, from the eigenvalue flows we see that there is in fact 
one more pair of would-be zero- modes of D st , also with 
positive chirality, giving index +4 in accordance with the 
Index Theorem 

A crucial property that would-be zero-modes of a lat- 
tice fermion formulation should have is robustness: they 
should not disappear under small deformations of the 
gauge held. This is assured in the present case if there is 
a clear separation between low- and high-lying crossing 
regions. Figure 3 shows the eigenvalue flow of H st (m) 
in the same background as Fig. 1 but over a larger m 
range. We see that the eigenvalue crossings only occur 
in a localized region around m = and in high-lying re- 
gions, \m\ > 9 in this case. The large separation between 
low-lying and high-lying crossings illustrates the robust- 
ness of the would-be zero-modes and index for staggered 
fermions. For comparison, Figure 4 shows the eigen- 
value flow of the hermitian Wilson operator H\y(m) in 
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FIG. 4: Wilson spectral flow in same background as Fig's 1,3. 



the same gauge field background. It has one low-lying 
positive-slope crossing in accordance with the Index The- 
orem, and the high-lying crossings are localized around 
to = 2 and m = 4 as expected on theoretical grounds 

am. 



As in the Wilson case [15[, separation between low- 
and high-lying crossing regions can be proved analyti- 
cally when the plaquette variables of the lattice gauge 
field satisfy the approximate smoothness condition 1 1 1 — 
U^ v {x)\\ < e for sufficiently small e [l7|: We derive in 18 1 
a bound of the form 



H s t{mf 



> 



m - 

1 - 



— K(m, d)e 
K(m, d)e 



for \m\ < 1 
for Iml > 1 



(8) 



(which is found to be saturated in the free field case 
where e = 0). The precise form of K(m,d) > is not 
important here, only the fact that it depends continu- 
ously on to, which ensures that Kq defined in the fol- 
lowing is finite. For any 61,62 with < 61 < 1 < 62 

b 2 

set Kq = m&x{K(m,d) \ \m\ < 62} and erj = j^, then 
the bound JHJ) implies H st {m) 2 > for b\ < \m\ < 62 
when e < eo in the plaquette condition. This shows that 
the separation between low- lying (|to| < 61) and high- 
lying (\m\ > 62) eigenvalue crossing regions can be made 
arbitrarily large by taking e > to be sufficiently small. 

The performance of the staggered index compared to 
the Wilson index was investigated in the numerical study 
with Q = l. For each randomly generated {r^x)} we ex- 
amined the would-be zero-modes and index as the rough- 
ness parameter S is increased. In both the staggered and 
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Wilson cases they were found to remain identifiable up 
to roughness levels 0.33 < 6 < 0.41 (which easily exceeds 
the limit S sa 0.25 at which the field-theoretic approach to 
the staggered fermion index was found to break down in 
Q). Furthermore, in all the backgrounds considered, the 
breakdown value of S for staggered was found to be the 
same as for Wilson up to differences A5 = ±0.02 (which 
favored staggered as often as Wilson). This suggests that 
the staggered and Wilson indexes are accessing the same 
topological content of the lattice gauge field. However, 
the computational cost was roughly twice as much in the 
Wilson case. This is as expected since d — 2 Wilson 
fermions have two spinor components. 

Finally we show that, analogously to the Wilson case, 
the staggered fermion index can be obtained as the index 
of the exact zero-modes of an overlap Dirac operator D ov . 
This is of practical as well as theoretical interest since the 
Wilson fermion index is usually calculated in practice as 
index(_D „). The role of 75 in the overlap construction 
is not to be played by T 5 here since it violates the 
required property 75 = 1 by an 0(a) term. Instead we 



(9) 



which corresponds to 75 ® 75 in the spin-flavor interpre- 
tation [lj| . We define the overlap Dirac operator with 
staggered fermion kernel by 



Dn 



H st (m Q ) 
y/H st (mo)' 



(10) 



As in the Wilson case [9j D ov has exact zero-modes, with 
definite chirality with respect to r 55 , and satisfying an 
index formula 



index (D ov ) 



H st (m ) 
y/H s t(m )' 



(11) 



which also follows from the general index formula in [16j] 
after noting that (fTO)) satisfies the Ginsparg- Wilson rela- 
16[ with 75 — > As in the Wilson case we take 



tion 



mo to be in the region in between the positive low-lying 
and high- lying eigenvalue crossings of H st (m) , then (| 1 1|) 
gives 



index(D ov ) — h'mdex(D s t) 



(12) 



To see this, note that -ff st (m)r55 = — T^^Hsti— m), which 
implies a symmetry in the eigenvalue flow: H st (m) and 
—H s t(—m) have the same spectrum (as seen in Figures 
1,2,3). It follows that H st (0) has symmetric spectrum; 
therefore (fTTj) is minus the spectral flow from m = to 
m = mo, which in turn is minus half the spectral flow 
from —mo to mo, i.e. half of index(_D s t). 

The staggered overlap operator introduced here is of 
independent interest as a new fermion formulation for 
Lattice QCD, and will be studied in a separate paper [lij . 



It has the remarkable feature of reducing the number of 
staggered fermion tastes by half (as reflected in the factor 
1/2 in (fl2")) L The physical fields turn out to correspond 
to the two continuum tastes with positive flavor-chirality 
under 1 (8 75, so that 75 <g> 75 chirality is the same as the 
physical 75 ® 1 for them. A new staggered version of 
domain wall fermions is also obtained |18| . 

In summary, staggered lattice fermions do maintain 
the important Index Theorem connection between gauge 
field topology and fcrmionic zero-modes, but in a way 
that was not realized previously. In the present d = 2 
study it was seen to perform as well as the Wilson in- 
dex, but with greater numerical efficiency. Future work 
should investigate the performance of the index for im- 
proved staggered fermions versus Wilson index in back- 
grounds generated in current Lattice QCD simulations in 
4 dimensions. 
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